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ABSTRACT 

Context. High resolution X-ray spectroscopy has revealed soft X-rays from high density plasma in Classical T-Tauri stars (CTTSs), 
probably arising from the accretion shock region. However, the mass accretion rates derived from the X-ray observations are consis- 
tently lower than those derived from UV/optical/NIR studies. 

Aims. We aim to test the hypothesis that the high density soft X-ray emission is from accretion by analysing optical accretion tracers 
from an X-ray selected sample of CTTSs in a homogeneous manner. 

Methods. We analyse optical spectra of the X-ray selected sample of CTTSs and calculate the accretion rates based on measuring 
Ha, H/3, Hy, He n 4686A , He i 5016A , He i 5876A , O i 6300A and He i 6678A equivalent widths. In addition, we also calculate the 
accretion rates based on the full width at 10% maximum of the Ho- line. The different optical tracers of accretion are compared and 
discussed. The derived accretion rates are then compared to the accretion rates derived from the X-ray spectroscopy. 
Results. We find that, for each CTTS in our sample, the different optical tracers predict mass accretion rates that agree within the 
errors, albeit with a spread of « 1 order of magnitude. Typically, mass accretion rates derived from Ha and He 1 5876A are larger than 
those derived from H/3, Hy and Oi. In addition, the Ha full width at 10%, whilst a good indicator of accretion, may not accurately 
measure the mass accretion rate. When comparisons of the optical mass accretion rates are made to the X-ray derived mass accretion 
rates, we find that: a) the latter are always lower (but by varying amounts); b) the latter range within a factor of » 2 around 2 x 10 -10 
M G yr~', despite the fact that the former span a range of « 3 orders of magnitude. We suggest that the systematic underestimation of 
the X-ray derived mass accretion rates could depend on the density distribution inside the accretion streams, where the densest part of 
the stream is not visible in the X-ray band because of the absorption by the stellar atmosphere. We also suggest that a non-negligible 
optical depth of X-ray emission lines produced by post-shock accreting plasma may explain the almost constant mass accretion rates 
derived in X-rays if the effect is larger in stars with larger optical mass accretion rates. 

Key words, accretion, accretion disks - circumstellar matter - stars: pre-main sequence - techniques: spectroscopic 

X i 1- Introduction in young stars and Br own D warfs (e.g. Natta et al. I2004I l2006l 
H Mohanty et al. I2003I I2005I I. Such studies are based on the hy- 
, The process of stellar mass accretion is an important aspect of pothesis that the UV veiling arises from the hot shock region 
star formation that is still to be fully understood. Not only is (e g Calvet & Gullbring [1998J, whereas the broad permitted 
accretion responsible for building up the young star to its final emission lines exhibited by CTTSs arise from the Mailing mag- 
mass, but it is also responsible for powering the mass outflows net ospheric flow (Muzerolle et al. IT9981 The Can lines have 
observed from such systems, which in turn remove the excess been found to be particularly good infa n tracers (Muzerolle et 
angular momentum, and prevent the star from spm-up. In ad- al rjgggj Mohanty et al. I20051 with the flux strongly corre- 
dition, understanding the mass accretion rate will significantly lated tQ me accretion rate> though other lines such as Ha and 

impact on the understanding of the inner disk, disk evolution u c „ 7 , s , , , , .• . 

\ , , . & . , _ ,, ' . , He i 5876A also show good correlations to the accretion rate 
and the eventual formation of planets. Stellar accretion can be /ZJ e „ „ , j fonnot 
, . „ • , „ ■ , j- , (Herczeg & Hillenbrand 2008). 
briefly summarised as follows: material passes from the enve- 
lope through the accretion disk - which is truncated at a radius For the last few decades it has been known that pre-main se- 
R m due to the strong magnetic field of the central star. At this in- quence stars have strong X-ray emission (Feigelson & DeCampli 
ner region of the disk, the material then flows along the star-disk |1981) . This X-ray emission was thought to have the same origin 
magnetic field lines (flux tubes) - at -free-fall velocity - onto as that of main sequence stars: low-density plasma (n e ~ 10 10 
the central star, where a strong shock is formed as the accreting cm 3 enclosed in coronal loop structures and heated to t emper - 
material impacts the stellar surface. atures of T ~ 10 6 - 10 7 K (Feigelson & Montmerle 1 19991 1. 

UV veiling, Ha, Pa/3, Bry, O i, Ca n and He i emission lines More recently, high resolution X-ray spectroscopy has revealed 

have all previously been used to probe the accretion scenario the presence of soft (E < 0.7 keV) X-ray emission, originat- 
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ing from high density (n > 10 11 crrT 3 ) plasma at temperatures 
of ~3 MK (see Telleschi et al. 120071 and references therein). It 
has been proposed that this soft X-ray emission is due to mass 
accretion. This interpretation is based on a simple model: as- 
suming the accretion flow has free-fall velocity v ~ 500 kms -1 , 
it becomes heated up by the shock (due to the impact of the 
accreting material with the stellar surface) at a temperature of 
T ~ (3 / 16)(/.miH\ 2 ) ~ 3 x 10 6 K, and then cools down radia- 
tively (Gullbri ng[T9M Calvet & Gullb ring[T998l Lamzin [T998l 
Gunther et al. 2007, Brickhouse et al. 120101) producing strong 
X-ray emission. 

Supporting this interpretation is the fact that this plasma 
component i) has never been observed in non-accreting stars, 
and ii) is too dense to have a coronal origin. This interpretation is 
also supported by time-dependent models of radiative accretion 
shocks in CT TSs (K oldoba et al. l2008l Sacco et al. [20081 [20101 
Orlando et al. 2010). In particular, Sacco et al. (2008) carried out 
a detailed hydrodynamic modeling of the interaction between 
the accretion flow and the stellar chromosphere, synthesizing the 
high resolution X-ray spectrum, as it would be observed with the 
Reflection Grating Spectrometers (RGS) on board the XMM- 
Newton satellite. They found an excellent agreement between 
predicted and observed X-ray spectra, supporting once again the 
idea that this X-ray emission originates mainly from accretion 
shocks. 

However, there are observational results that are difficult to 
reconcile with this framework. In particular, the mass accretion 
rates derived from X-rays are usually underestimated if com- 
pared to accretion rates derived from optical and UV data. As an 
example, for TW Hya, X-rays indicate M ~ 1 x 10 _11 M o yr _I 
(Stelzer et al. 2004), while Ha and UV provide 5 x 10 10 M yr" 1 
(Muzerolle et al. 2000). The same mismatch holds for BP Tau: 
M ~ 9 x 10 -10 Moyr 1 from X-rays (Sch mitt et al. l2005l and 
3 x 10~ 8 M yr" 1 from UV (Gullbring et al. [1998b . To date, only 
3 such comparisons have been made in the literature (TW Hya; 
BP Tau; MP Mus, Argiroffi et al 2009), moreover, each of these 
comparisons has used differing UV/optical accretion measures 
and differing means of calculating the X-ray accretion rate. In 
order to confirm that the soft X-ray emission arises from the ac- 
cretion process, and understand the discrepancies between the 
optical and X-ray accretion rates, we have homogeneously anal- 
ysed optical data for all the CTTSs for which high-resolution 
X-ray spectra, with good S/N in the [0.5, 1.0] keV range, have 
been gathered. We compare the X-ray derived mass accretion 
rates to those measured from the HB, Hei 5876A, Oi 6300A, 
Ha equivalent width and Ha full- width at 10%. 

The structure of the paper is as follows. In section 2 we 
present our sample of CTTSs previously observed with Chandra 
and XMM. We briefly discuss the parameters of each CTTS 
(summarised in Table[TJ, along with the previously published X- 
ray data. In section 3 we discuss the data reduction and analysis 
of the optical data (subsection 3.1), and the derivation of mass 
accretion rates from the X-ray data (subsection 3.2). In section 
4.1 we draw comparisons between the different optical tracers 
of accretion and discuss the mass accretion rates derived from 
these tracers. In section 4.2 we analyse and discuss the variabil- 
ity of TW Hya. In section 4.3 we compare the optical and X-ray 
mass accretion rates and discuss possible scenarios to explain 
our results, and in section 5 we present our conclusions. 

2. The Sample 

The sample consists of all the CTTSs currently observed with 
high-resolution X-ray spectroscopy (either with Chandra, or 



RGS on XMM-Newton) and for which O vn triplet lines have 
been measured, and which also have high resolution optical 
echelle spectroscopy available in the data archives. The X-ray 
data have previously been published, with the published fluxes of 
the O vn used to derive the mass accretion rate for each star. The 
parameters for each source, used in our calculations are given in 
TableQ] Below are brief descriptions of each source. 

Hen 3-600 is a multiple star system and a member of the 
TW Hydra Association, at a distance of 45 pc. This is an aver- 
age distance taken from Huenemoerder et al. (120071) . who used 
the photometric distance from Kastner et al. ( 1997 ) and the typi- 
cal distance adopted for stars belonging to the TW Hya associa- 
tion. The primary components (A and B) of this system are sep- 
arated by l'/4, and component A is surrounded by a dusty disc 
(Jayawardhana et al. 119991) . Component A is also inferred (by 
Huenemoerder et al. 120071) to be almost pole-on, from a large 
mid-IR excess and a negligible optical reddening (based on a 
B-V = 1.52 which is nearly that of an unreddened M3 photo- 
sphere; Johnson [1966I) . Huenemoerder et al. (120071 ) discuss the 
X-ray Chandra observations of Hen 3-600, and find high density 
plasma, and a larger 'soft excess' for component A, than B. 

TW Hya is one of the closest known CTTSs, at a distance 
of only ~ 56 pc (Wichmann et al. I1998I ). It has a mass of 0.7 
M G , a radius of 1 R (Muzerolle et al. 2000) and is orientated so 
that it is seen nearly pole-on (Kastner et al. 1997). Previous X- 
ray data have been published by Kastner et al. (2002j), Raassen 
(2009), Brickhouse et al. d20T0l ) and Stelzer & Schmitt d2004t , 
who found a strong soft X-ray emission produced by high den- 
sity plasma. 

RU hup is located at a distance of 140 pc (Hughes et al. 
1993). It has a mass and radius of 0.8 M and 1.7 R respec- 
tively (Stempels & Piskunov 120021 ). it suffers little absorption 
(Herczeg et al. [2005t and is suggested to be viewed almost pole- 
on (Stempels & Piskunov 120021 ). XMM data are discussed by 
Robrade & Schmitt (2007), who find cool, high density plasma 
which they conclude indicates an accretion shock origin. 

BP Tau is a CTTS in the Taurus-Auriga molecular cloud at a 
distance of 140 pc (Kenyon et al. 2008). It has a spectral type of 
K 7, and a mass of 0.8 M G (Kenyon & Hartmann "l995l Giidel et 
al. 120071 ). The XMM-Newton spectra were previously discussed 
by Schmitt et al. d20"05i ) and Robrade & Schmitt d2U06l ). They 
report the presence of high density soft X-ray emission from the 
O vn triplet, suggesting this could be due to accretion shocks. 

V4046 Sgr is a nearby spectroscopic binary CTTS that is 
isolated from any dark cloud or molecular cloud, it has negli- 
gible extinction, and there is evidence for a circumstellar disk 
(Hutchinson et al. 1990). The period is well determined to be 
2.4213459 days (Stempels & Gahm 12004b and the separation is 
estimated to be ~ 10 R (Quast et al.|2000). The observed spec- 
tral energy IR distribution is consistent with the disk having an 
inner radius of 1.8AU, therefore the disk is circumbinary. The 
inclination of the system is ~35^45°. Chandra data of V4046 
Sgr is presented by Gunther et al. (2006 ), who found the X-ray 
emission is due to high density plasma. 

MP Mus is a Kl IVe type star located in the Lower Centaurus 
Crux (LCC) association at a distance of ~ 86 pc. It is known to 
have a dusty disk (Mamajek et al. 2002; Silverstone et al. 2006) 
with a dust mass of ~ 5 X 10~ 5 M (Carpenter et al. 120051 ) . 
Batalha et al. ( 119981 ) studied the variability in B, V, R and I 
bands, and found that it has a low variability for a CTTS. XMM 
data for MP Mus have previously been published by Argiroffi 
et al. ( 12007b . who found evidence for a high density plasma re- 
sponsible for the soft X-ray emission. 
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Table 1. List of sources and their adopted parameters. 



Name 


RA 


Dec 


Dist. 


Mass 


Radius 


Age 


Sp Type 


Mag. 


A v 


Incl. 


Binarity" 


References 




(J2000) 


(J2000) 


(pc) 


(Mo) 


(Rg) 


Myr 




(Johnson V) 


(mag) 


C) 






Hen 3-600 


11 10 28.1 


-37 31 51 


45 


0.2 


0.9 


8 


M4Ve/M4Ve 


12.04 


0.7 


~0 


BS 


1, 2, 3, 4, 5, 6 


TWHya 


110151.9 


-34 42 17 


56 


0.7 


1.0 


8 


K7 


11.27 





7 


ss 


1, 6, 7, 8, 9, 10, 11, 12 


RULup 


15 56 42.3 


-37 49 15 


140 


0.8 


1.7 


0-3 


K7 


11.55 


0.1 


24 


ss 


7, 9, 13, 14, 15, 16 


BP Tau 


04 19 15.8 


+29 06 27 


140 


0.8 


2.0 


1.9 


K7 


12.13 


0.5 


45 


ss 


9, 17, 18, 19, 20, 21, 22 


V4046 Sgr 


18 14 10.5 


-32 47 34 


72 


0.86 


1.16 


12 


K5Ve/K7Ve 


10.69 





35 


BS 


9, 23, 24, 25 


MP Mus 


13 22 07.5 


-69 38 12 


86 


1.2 


1.3 


6-7 


KUVe 


10.44 


0.17 


32 


SS 


9, 26, 27, 28 


V2129 Oph 


16 27 40.3 


-24 22 03 


120 


1.35 


2.4 


2 


K5 


12.28 


0.3 


45 


ss 


3, 9, 29, 30 


T Tau N 


04 21 59.4 


+ 19 32 06 


140 


2.4 


3.6 


2.7 


K0 


9.88 


1 


13 


TS 


8, 15, 16, 20,31 



SS = single star; BS = binary system; TS = triple system. 



Refer ences. (1) Muzero lle et al. 125551 ( 2) Tor res et al. 125551 (3) Geoffray & Monin l255Tl ( 4) Huenemoerder et al .l2007l (5) Kastner et al. [l997l (6) So ng et al. 120031 (6) Herczeg & Hillenb rand 
120081 (7) Qi et al. l2004l (8) Lasker et al. 120081 ( 9) Wic hman et al. [T998l (10) R aassen 120091 (1 1) Brickho use et al. l20T0l ( 12) Ro brade & Schmitt 12507] (13) Stempels, Gah m & P etrov l2007l (14) 
Stempels & Pisk unov [2002 (15) Comer on et a l. 120081 (15) Giidel et al. 120071 (16) Grankin et al. 120071 (17) Costa et al. 120001 (18) Donati et al. 2008; (19) Kenyon et al. 20081 (20) Kenyon & 
Hartm ann [T995l (21) St empels & Gahm l2004l ( 22) Hu tchinson et al. [19901 (23) Kastner et al. l2008l (26) Mamajek et al. l2002l (27) Cortes et al. |2009l (28) Torres et al.[2008; (29) Donati et al. 
120071 (30) Wilking et al. 125531 (31) Bertout et al. [T^88l 
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Table 2. Optical archival data information for our sample of CTTSs. 



Name 


Telescope 


Instrument 


Observation Date 
(yyyymmdd) 


Exposure time 

(s) 


Hen 3-600 


ESO 2.2m 


FEROS 


20040512 


900 


TWHya 


ESO 2.2m 


FEROS 


20070426 


900 


RULup 


VLT 


UVES 


20050814 


900 


BP Tau 


TNG 


SARG 


20071220 


3600 


V4046 Sgr 


VLT 


UVES 


20050815 


900 


MPMus 


ESO 2.2m 


FEROS 


20060418 


3000 


V2129 Oph 


VLT 


UVES 


20020417 


900 


T Tau 


TNG 


SARG 


20061130 


3600 



V2129 Oph is the brightest CTTS in the p Oph cloud, at a 
distance of 120 pc (Lombardi, Lada & Alves (2008 )). It is a dis- 
tant binary, with its very low-mass companion (possibly a brown 
dwarf) being about 50 times fainter than the CTTS in the V band 
(Geoffray & Monin 2001). The CTTS is inferred to have a mass 
of 1.35 M and a radius of 2.4 R . Chandra data have recently 
been obtained (Argiroffi et al. in prep) but not published yet. We 
include the optical analysis in this paper for easy comparison 
once the X-ray data is published. 

T Tau N is another CTTS located in Taurus-Auriga. It is op- 
tically visible part of a triple system, along with T Tau Sa and 
Sb which are the "infrared companions". It has a spectral type of 
K (Ken yon & Hartmann [l9951 ) and a mass of 2.7 M (Giidel et 
al. I20071 I. It is oriented such that it is seen almost pole-on (Solf 
& Bohm [1999l Akeson et al. [19981 . X-ray d ata fro m XMM and 
Chandra have been discussed in Giidel et al. (2007) who found a 
'soft excess' but no evidence of high density plasma. 



the observed V magnitude we take into account the continuum 
excess due to the accretion process. 

The mass accretion rates were calculated following the em- 
pirical relation between accretion and line luminosities found 
and described by Herczeg & Hillenbrand (2008): 



log L acc = a + b log L a 



(1) 



where the coefficients a and b (listed in Tab. [4]) have been cal- 
culated - by Herczeg & Hillenbrand ( 2008) - by comparing op- 
tical emission line fluxes and accretion luminosities measured 
from the UV continuum excess. Errors on the coefficients a and 
b depend on the scatter of the line luminosity values around the 
best fit relation (from a large sample of sources), and therefore 
take into account several effects, such as optical depth and any 
contribution to the line emission from stellar outflows. 

Once the accretion luminosity, L. dcc , has been calculated, it is 
possible to calculate the mass accretion rate via: 



3. Data Analysis 

3.1. Optical Data Analysis 

We retrieved high resolution echelle spectra spanning the optical 
wavelength range from the archives. No single instrument had 
observed the entire sample (in part, due to their locations). Data 
from three instruments, FEROS, on the 2.2 m ESO telescope in 
Chile, UVES, on the VLT in Chile and SARG on the TNG in La 
Palma were used in this analysis. Details of the telescope, instru- 
ment, observation date and exposure times for each source can 
be found in Table [2] Relevant calibration files were downloaded 
with each dataset and used in the data reduction process. 

The data from the three instruments were reduced in the stan- 
dard manner using IRAF routines to de-bias, flatfield and extract 
the spectra. Wavelength calibrations were carried out based on 
ThArNe arcs for FEROS, ThAr arcs for UVES and Thorium arcs 
for SARG. 

Line luminosities, Lij ne , have been calculated from the equiv- 
alent width of each emission line (see table 0. The contin- 
uum fluxes have been estimated from the V magnitudes of the 
stars, following a method similar to that used by Mohanty et 
al. (2005) and Dahm (2008). Specifically, the continuum flux at 
5500 A was calculated using the extinction corrected observed 
V magnitude, the distance, and the photometric zero point of 
the Johnson V filter (at the central wavelength 5500 A) of 3.75x 
10~ 9 erg cm -2 s _1 A -1 (Mitchell & Johnson [19691 ) . The contin- 
uum flux at the wavelength of each of the analysed emission 
lines also required correction for the stellar spectral response - 
which was measured using the Pickles Spectrophotometric Atlas 
of Standard Stellar Spectra (Pickles 119851 ) in conjunction with 
the observed stellar spectral type. It is worth noting that using 



M= 1 



R* 

Rtn 



R* 
GM, 



1.251* 



R* 
'GMl 



(2) 



where (1 -R*/Ri n ) 1 w 1.25 is estimated by assuming the accret- 
ing gas falls onto the star from the truncation radius of the disk, 
R in * 5R* (Gullbring et al. [19981 . 

We also calculated the mass accretion rate based on the Ha 
full width at 10%, using the following equation (Natta et al. 
2004): 



log M * -12.9(+0.3) + 9.7(+0.7) x 10~ 3 Hal0% 



(3) 



where Hffl0% is the Ha 10% full width in kms 1 and M 
is in M yr -1 . The derived mass accretion rates, along with a 
(weighted) mean optical accretion rate, (Mo p t), are listed in 
Table[5] The mean optical accretion rate was calculated using all 
the mass accretion rates calculated via Eqs.Q] and [2] and there- 
fore does not include the mass accretion rate derived from the 
Ha 10% width (see section l4~Tl for more details). 

3.2. X-ray Data Analysis 

It is also possible to infer the mass accretion rate M of CTTSs 
from their soft X-ray emission (e.g. Schmitt et al. [2005). This 
method for deriving M is based on the assumption that the whole 
soft X-ray emission is due to accreting material and not to coro- 
nal plasma, and that the post-shock zone can be described as 
a stationary isothermal slab of plasma at constant velocity and 
density. Recent hydrodynamical simulations (Sacco et al. 2008; 
120101 ) have demonstrated that, even if the post-shock zone is 
not stationary (and characterized by quasi-periodic shock oscil- 
lations), the time-averaged properties of the post-shock zone are, 
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Table 3. Measured equivalent widths for the accretion tracers" , along with the measured full width at 10% for the Ha line. 



Source 


Hy 4340A 


Hen4686A 


H6 4861A 


Hei5016A 


Hei5876A 


Oi6300A 


Ha 6563A 


Ha FW10% 


Hei6678A 




A 


A 


A 


A 


A 


A 


A 


kms -1 


A 


Hen 3-600 


-7.51 




-7.93 


-0.32 


-0.96 


-0.37 


-21.81 


282.29 


-0.37 


TW Hya 


-18.49 


-0.48 


-31.9 


-0.32 


-3.16 


-0.54 


-141.58 


407.55 


-0.7 


RULup 


NA 


NA 


-28.34 


-9.25 


-4.55 


-1.17 


-68.95 


578.42 


-1.88 


BPTau 


NA 


-0.6 


-19.82 


-0.27 


-1.08 


-0.37 


-58.63 


413.41 


-0.31 


V4046 Sgr 


NA 


NA 


-5.97 




-0.41 


-0.12 


-32.66 


515.19 




MP Mus 


-1.34 




-3.16 




-0.3 


-0.23 


-23.47 


502.9 


0.16 


V2129 Oph 


NA 


NA 


-2.54 




-0.33 


-0.13 


-12.07 


274.24 


-0.04 


TTau 


NA 


-0.1 


-12.05 


-0.6 


-0.55 


-0.7 


-42.45 


419.71 


0.08 



Negative values represent emission. Note some sources show lines in absorption. There are no measurements where the spectra show no 
indication of a line present. NA indicates where the wavelength of the line was not covered in the observation. 



Table 4. Coefficients used in Equation[T| for each emission line, 
from Herczeg & Hillenbrand (l2U08l l 



Line 


Wavelength 

(A) 




a 


b 


Hy 


4349 


3.0 


±0.2 


1.24 ±0.04 


Hen 


4686 


3.7 


±0.5 


1.04 ±0.08 


H8 


4861 


2.6 


±0.2 


1.22 ±0.05 


He i 


5016 


3.3 


±0.3 


1.04 ±0.05 


He i 


5876 


5.3 


± 0.7 


1.46 ±0.12 


Oi 


6300 


2.8 


±0.8 


0.96 ±0.16 


Ha 


6563 


2.0 


±0.4 


1.2 ± 0.11 


He i 


6678 


5.2 


±0.8 


1.37 ±0.13 



in general, well-described by the stationary model. Sacco et al. 
(120101 1 have also shown that the higher the post-shock tempera- 
ture the larger the discrepancies between the stationary and the 
hydrodynamical models. This is mainly due to the fact that the 
thermal conduction (which is not taken into account in the sta- 
tionary model and that is more efficient for higher post-shock 
temperatures) drains energy from the shock-heated plasma to 
the chromosphere through a thin transition region, and acts as 
an additional cooling mechanism (see also Orlando et al. 2010). 

Letting ra an d v o represent the density and velocity of the 
pre-shock material, and n\ and Vi the corresponding post-shock 
values, then pre- and post-shock quantities are linked by these 
relations: 



«i = 4«o 



V! = ~V0, 



r, = 



3 yUffl H 2 

i6ir Vo 



(4) 



where T\ is the post-shock temperature, fi is the average mass 
per free particle (/i « 0.5), and we assume a strong shock regime. 
The volume occupied by the hot post-shock plasma is defined by 
the stream cross-section, A, and by the distance, I, covered by the 
plasma itself before cooling down. Hence, when r is the cooling 
time, / is given by Vit, and the emission measure - EM - of 
the shock-heated plasma is n 2 x A\\T. Therefore it is possible to 
derive the stream cross-section A, and hence the mass accretion 
rate, from the EM value measured from the soft X-ray emission 
via: 



. EM EMP{T\) 

M = juranniv^ = jura H = yum H — r— 

n\T ikl i 



(5) 



where r = 3kT l{n\P{T)), with P(T) indicating the plasma radia- 
tive losses per EM unit. We used the above formula to obtain M 
from the X-ray data. In particular for each CTTS: 1) we evaluate 
vq, and hence vi, assuming a free fall velocity from a distance of 



5/?* (to correspond with the optical analysis); 2) we evaluate the 
post shock temperature T\ from vo (see Eq.|4j>; 3) the knowledge 
of T\ allows us to derive the plasma EM from the observed flux 
of the O vii resonance line located at 21.60 A via: 



EM : 



^ovu 



G(r,)A c 



(6) 



where G(T[) is the emissivity function of the O vn resonance line 
at 21.60 A, and Ao is the oxygen abundance. 

We chose to rely our M estimations on the flux of the O vn 
resonance line because it is mostly produced by plasma at 2 MK, 
which is the typical expected temperature for the plasma heated 
in the accretion shock, and because its flux measurement is avail- 
able for almost all CTTSs observed with high resolution X-ray 
spectroscopy. 

Note that this method allows us to derive M independently of 
the knowledge of the plasma density n \ . On the other hand this 
method depends on the plasma metallicity, which is needed to 
calculate the radiative losses P(T\) and to link the O vn line flux 
to the corresponding EM. It is worth noting that X-ray obser- 
vations provide us calibrated spectra, making the line flux mea- 
surements a simple and solid procedure, especially for the O vn 
resonance line at 21.6 A that is an isolated line, and hence its 
flux estimation is not affected by problems due to line blending 
or wrong continuum placement. We decided, therefore, to infer 
the M values, with the above method, starting from published 
values of the fluxes of the O vn resonance line. 

In Table|6]we list, for each CTTS, the derived accretion rates 
together with the main parameters used. 

4. Results & Discussion 

4.1. Comparison of optical derived accretion rates 

Table [5] shows the calculated mass accretion rates for each star 
based on the different emission lines analysed. Although the 
sources were observed at different epochs, all the emission lines 
were observed simultaneously for each star. Therefore one might 
naively expect the different tracers to lead to very similar mass 
accretion rates (if each line is tracing the same accretion phe- 
nomena). The top panel of Fig.[T]shows the different mass accre- 
tion rates calculated for each star in the sample (in order of in- 
creasing stellar mass). We find that all the different tracers agree 
within the errors, albeit with a spread of typically ~1 order of 
magnitude. We see no relation between mass accretion rate and 
the mass of the stars. Also, it should be noted, that there is no 
apparent relation between the mass accretion rate and the age of 
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Table 5. Derived optical mass accretion rates based on the equivalent widths and Ha full width at 10% measurements listed in Table [3] Also listed are the optical weighted mean 
accretion rates, (Mopt), based on all the accretion rates except Hal0% width. 



Source 


logM ? 

c Hy 4340A 


logM j 

c Hen 4686 A 


logM s 

° H/3 486lA 


logM j 

° He I 5016A 


logM i 

c He I 5876 A 


logM % 

B Ha 6563A 


log M Ka io% 


logM ? 

c O I 6300A 


loaM s 

c He I 6678 A 


log<M 0pt > 




Meyr- 1 


M e yr- 1 


M s yr' x 


M Q yr- 1 


Meyr- 1 


M s yr' x 


M Q yr~ l 


M Q yr l 


M Q yr l 


M Q yr-' 


Hen 3-600 


-9.75 ±0.28 




-9.61 ±0.30 


-9.52 ± 0.42 


-9.02 ± 0.94 


-8.89 ±0.55 


-10.15 ±0.61 


-9.27 ± 1.19 


-8.89 ± 1.06 


-9.53 ±0.17 


TWHya 


-9.31+0.26 


-9.17 ±0.66 


-9.05 ± 0.27 


-9.82 ± 0.41 


-8.38 ±0.88 


-8.46 ±0.50 


-8.94 ± 0.78 


-9.39 ± 1.15 


-8.91 ± 1.03 


-9.17 ±0.15 


RULup 






-8.06 ±0.25 


-7.38 ± 0.35 


-6.91 ± 0.82 


-7.79 ±0.47 


-7.28 ± 1.06 


-8.20 ± 1.03 


-7. 16 ±0.94 


-7.77 ±0.18 


BPTau 




-8.14 ±0.62 


-8.26 ±0.25 


-8.97 ± 0.39 


-7.86 ± 0.86 


-7.89 ± 0.48 


-8.88 ± 0.79 


-8.68 ± 1.09 


-8.26 ± 1.01 


-8.34 ±0.17 


V4046 Sgr 






-9.34 ±0.28 




-9.07 ± 0.91 


-8.79 ±0.51 


-7.89 ± 0.95 


-9.62 ± 1.17 




-9.22 ±0.23 


MPMus 


-9.61+0.26 




-9.32 ± 0.27 




-8.89 ± 0.90 


-8.73 ±0.50 


-8.01 ± 0.92 


-9.18 ± 1.11 




-9.35 ±0.17 


V2129 Oph 






-9.76 ±0.29 




-9.19 ±0.92 


-9.28 ±0.55 


-10.23 ±0.60 


-9.56 ± 1.18 


-9.99 ± 1.11 


-9.64 ± 0.24 


TTau 




-7.90 ± 0.60 


-7.22 ± 0.23 


-7.50 ± 0.35 


-6.98 ±0.81 


-7.16 ±0.44 


-8.82 ± 0.80 


-7.68 ± 0.97 




-7.32 ±0.16 
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Fig. 1. Top Panel: Plot of the different mass accretion rates from different accretion tracer emission lines, for each of the stars in 
the sample. The stars are in order of increasing mass, showing no relation between mass and mass accretion rate. Bottom Panel: 
Plot showing the ratio of each mass accretion rate estimate to the optical mean of the accretion rate for each star. The different 
symbols/colors represent different accretion rate tracers (see upper left corner in the top panel). 



Table 6. Observed X-ray parameters and derived X-ray mass accretion rates . 



Name T O, Ne, Fe P(T) logLovn log£M logM References for L vn 

(MK) (l(T 23 cm 3 erg s~') (ergs -1 ) (cm -3 ) (M yr-1) and abundances 

Hen 3-600 TO 0.4, 1.2, 0.2 5/7 27.5 ± 0.2 52.7 ± 0.2 -9.53 ± 0.26 I 

TWHya 3.1 0.2,1.8,0.3 2.2 28.9 ±0.1 53.4 ±0.1 -9.74 ± 0.05 2 

RULup 2.1 0.6,1.4,1.1 7.2 28.9 ±0.1 52.8 ±0.1 -9.67 ±0.13 3 

BPTau 1.8 0.6,1.5,0.3 4.7 28.9 ±0.1 52.8 ±0.1 -9.78 ± 0.08 4,5 

V4046Sgr 3.3 0.1,1.0,0.1 1.6 28.5 ±0.1 53.5 ±0.1 -9.84 ±0.11 6,7 

MPMus 4.1 0.1,0.4,0.1 1.3 28.6±0.1 53.8 ± 0.1 -9.70±0.06 8 
V2129 0ph 

TTau 2.9 0.4,0.8,0.3 2.4 29.4 ±0.1 53.6 ±0.1 -9.47 ±0.11 9,10 



{a) Chandra data for V2129 Oph have recently been obtained (Argiroffi et al. in prep) but not published yet. 

Refere nces. ( 1) Huenemoerder et al. [20071 (2) Stelzer & Schmitt, l2004l (3) Robrade & Schmitt l2007l (4) Schmitt et al. [20051 (5) Robrade & 
Sch mitt 120061 (6) Gunther et al. l2006l (7) Argiroffi C, private communication; (8) Argiroffi et al. [2009l (9) Giidel & Telleschi [20071 (10) Giidel 
etal. [2007l 
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Fig. 2. The observed Ha line profiles for our sample. The data show many different line morphologies, with some very asymmetric 
lines e.g. V2129 0ph. 



the stars (see table [TJ - for example, V2129 Oph has an age of 
2 Myrs, whereas T Tau has an age of 2.7 Myrs and their mass 
accretion rates are vastly different. V4046 Sgr has an age of 12 
Myrs and MP Mus has an age of 6-7 Myrs, and yet they have 
similar mass accretion rates. The bottom panel of Fig. [TJ shows 
the distribution of the difference between each mass accretion 
rate estimate and <Mo p t), the optical mean mass accretion rate. 
In general, Ha flux and Hei 5876A lead to higher mass accre- 
tion rates than H/3, Hy and O i. 

Ha 10% mass accretion rates are often at the extremes of 
the range of estimates, but this is not surprising given the Ha 
line profiles are dependent on inclination, rotation effects (see 
Mohanty et al. 2005 for details), winds, and outflows. Fig. [2] 
shows the Ha line profiles for our sample, some of which are 
very asymmetric. The Ha emission from Hen 3-600 and MP 
Mus both have narrow peaks, but with asymmetric broad wing 
emission. TW Hya, RU Lup and BP Tau all have broader peaked 
emission, albeit with some evidence for self-absorption. V4046 
Sgr has an Ha line profile which has a broad peak, with very 
broad wing emission. V2129 Oph and T Tau both show clear 
evidence for absorption in the Ha emission, with V2129 Oph 
having a very asymmetric profile with strong absorption of the 
blue-shifted wing/no blueshifted component. Interestingly, we 
see that the accretion rates derived from Ha 10% emission is 
to the lower part of the range of mass accretion rates for Hen 



3-600, BP Tau, V2129 Oph and T Tau, at the upper end of the 
range for RU Lup, V4046 Sgr and MP Mus and in the middle of 
the range for TW Hya (see top panel of Fig. [TJ. The asymmetry 
of the emission line profile from V2129 Oph would clearly lead 
to the under-estimation of the mass accretion rate for that object, 
but interestingly, the under-estimation is not extreme - the dif- 
ference between the Ha 10% mass accretion rate and the optical 
mean mass accretion rate being only -1.7 x 10~ 10 M yr~'. In 
the cases of V4046 Sgr, MP Mus and T Tau, these estimates are 
clear outliers from the rest of the accretion rate estimates (high 
estimates for the first two - as predicted from the line profiles, 
and lower for T Tau - which has a 'normal' line profile, albeit 
with some absorption, but this does not affect the 10% intensity 
level). Due to the fact that Ha 10% mass accretion rate estimates 
are outliers in almost half our sample, along with the fact that 
inclination, rotation, winds and outflows may affect these mass 
accretion rates, we have excluded the Ha 10% mass accretion 
rate estimates from the calculation of the optical mean mass ac- 
cretion rate listed in Table [5] 

4.2. Variability of TW Hya 

TW Hya is the only source within the sample to have been 
observed over a period of months (from March 2007 through 
to July 2007) on both short (nightly) and long (monthly) 
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Fig. 3. The variability of TW Hya. The different color/symbols 
represent different accretion rate tracers (see upper right corner). 
Note the X-axis timescale is not linear. The variability in accre- 
tion rate is observed to be small in comparison to the range of 
accretion rates calculated from the different tracers. 



timescales, using the same instrumental set-up. We have anal- 
ysed 14 nights of data over a 5 month period for the purpose 
of studying the accretion variability. Our results are shown in 
Fig. [3] We find the variability for a given tracer small compared 
to the ~ 1 order of magnitude spread in mass accretion rates 
from the different tracers. Oi, Ho- and Hei 5876A are the trac- 
ers with the least variability, whilst H/3 and Hy have the largest 
variability. It could be argued that this small variability is ex- 
pected given TW Hya is viewed almost pole-on, and as such it 
should be possible to observe all active accretion regions about 
the facing pole at all times, therefore eliminating large changes 
in mass accretion rates as the hotspots rotate out of view. The 
variability of mass accretion rates are not necessarily due only 
to the inclination angle of the star-disk system, but may be due 
to time-dependent accretion rates. The small amount of variabil- 
ity (~ 0.25 order of magnitude) we observe in our TW Hya data 
may be due to time-dependent accretion rate. 

Alencar & Batalha (2002) studied the TW Hya accretion rate 
for 1.5 yrs, and found the mass accretion rate to be variable be- 
tween 10~ 9 and 10~ 8 Mq yr~' using the Na d line profile. Whilst 
these results fall within the range of mass accretion rates we 
measured using different tracers, they suggest a larger variabil- 
ity. This discrepancy may be due to the method used to estimate 
the mass accretion rate or due to the longer timeperiod covered 
in their analysis. Rucinsky et al. (2008) analysed two datasets 
- from 2007 and 2008 - and found a photometric variability of 
AV < 0.5 mag on timescales of ~ 3.5 days in the 2007 data. 
This however was no longer observable in the 2008 data, which 
showed variability (again AV < 0.5 mag) on timescales of 2-9 
days, which they suggest may be due to the orbital decay of ac- 
creting gas. In addition, they find 'spikes' in the accretion rate 
lasting only a fraction of a day. They conclude that whilst accre- 
tion contributes to this variability, it is not the sole cause. 

Our optical and X-ray observations are not coeval, and so it 
is important to study the variability of these stars due to accre- 
tion to know if the their variability will limit our ability to com- 
pare mass accretion rates derived from non-coeval observations. 
We find TW Hya, a pole-on star, to have a variability within 
our errors. Furthermore, all the sources in this sample have in- 
clinations < 45°, with TW Hya, Hen 3-600, RU Lup and T Tau 
all having inclinations much closer to pole-on, meaning any ac- 
cretion hotspots would be visible all/most of the time, limiting 
variability due to observability of the hotspot. Whilst we cannot 



rule out time dependent mass accretion rates for the stars other 
than TW Hya in our sample, we measure a small variability in 
our TW Hya data, and even the larger variabilities suggested in 
the literature fall within our errors. 



4.3. Comparison of optical/x-ray derived accretion rates 

We compare the optical derived mass accretion rates (Mopt) to 
the X-ray derived mass accretion rates Mx-my in Fig. [4] The most 
striking result is that Mx-my for our sample ranges within a factor 
* 2 around 2x 10~ 10 Mp yr _I (see Fig.|]and Table©, despite the 
fact that the range of (Mopt) spans almost 3 orders of magnitude. 
In addition, as suggested by other studies in the literature, the X- 
ray mass accretion rates are always lower than the corresponding 
optical mass accretion rates with the discrepancy increasing for 
increasing (Mo pl ). 

The above results can be interpreted in the light of the find- 
ings of Sacco et al. (120101 1 who used a detailed hydrodynamic 
model to investigate the observability of accretion shocks in 
X-rays in a wide range of physical conditions of the accretion 
stream. They considered the absorption from the optically thick 
plasma of the stellar chromosphere on the X-ray emission arising 
from the shock-heated material. They found that the absorption 
strongly depends on the accretion stream properties and influ- 
ences the observability of the post-shock accreting plasma in the 
X-ray band. In particular, observable X-ray emission from ac- 
cretion shocks is expected in streams with densities in the range 
10" <n e < 10 12 cirT 3 and velocities in the range 400 < v < 600 
km s _1 . In fact, denser streams produce post-shock zones that are 
rather narrow and deeply rooted in the chromosphere, causing 
their X-ray emission to be strongly absorbed. Streams with ve- 
locities below 400 km s _1 produce shocked plasma with temper- 
atures around 1 MK or even lower and, therefore, hardly observ- 
able in X-rays. The obvious consequence is that, whilst accretion 
shocks are almost always observable in the optical band, only a 
small fraction of them can be revealed in X-rays, leading in gen- 
eral to an underestimate of the X-ray derived mass accretion rate 
and, therefore, to the observed discrepancy with the optical de- 
rived mass accretion rates (see also discussion in Orlando et al. 
l20T0l and Sacco et al. [20T0l . 

Note that the above argument is applicable even if only one 
stream is present on the star. In fact, for typical densities and ve- 
locities of accretion streams and stellar magnetic field strengths 
(as derived from observations of CTTSs; Johns-Krull 2007), in 
general the plasma (3 (i.e. the ratio gas pressure / magnetic pres- 
sure) is expected to be <K 1 (see Sacco et al. 2010). In these 
conditions, an accretion stream is a bundle of fibrils each inde- 
pendent of the others due to the strong magnetic field which pre- 
vents mass and energy exchange across magnetic field lines (e.g. 
Orlando et al. 120 101 and Sacco et al . 120 1 01 > . If the stream is struc- 
tured in density (as suggested, for instance, by 3D MHD models 
of the star-disk system carried out by Romanova et al. |2004), 
only fibrils with density and velocity in the range of observabil- 
ity defined by Sacco et al. (2010) are expected to produce de- 
tectable X-rays, unless they are located close to the centre of the 
stream cross-section where they suffer strong absorption from 
the accretion column itself. 

Interestingly, Hen 3-600 is the lowest mass star in our sample 
and has a close agreement between the optical and X-ray mass 
accretion rates, suggesting it has very few high density streams 
(or fibrils) - hardly detectable in X-rays - or only low density 
accretion streams (fibrils) - observable in X-rays - meaning al- 
most everything is seen in the X-ray, as well as in the optical. At 
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Fig. 4. Plots of the optically derived mass accretion rates versus the X-ray derived mass accretion rates. Each panel shows the optical 
accretion rate derived from a different tracer, indicated in the lower right corner of the plot. The CTTSs are colour coded (see lower 
left corner of each panel). We plot only those cases for which we have measured the mass accretion rate for all the stars in our 
sample (see table[5]). The bottom right panel shows a plot of the optical mean mass accretion rate plotted against the X-ray accretion 
rate. There is a very small range, a factor « 2, in X-ray accretion rate calculated for the entire sample, whereas the optical mass 
accretion rates span a range of ~ 3 orders of magnitude. 



the other end of the scale, RU Lup could be interpreted to have 
a much larger number of high density streams (fibrils). 

If, on one hand, the discrepancy between optical and X- 
ray derived mass accretion rates may be easily explained in the 
framework of accretion shock models as discussed above, on the 
other hand, more puzzling is the result that Mx-ray is almost the 
same for all the stars of our sample. This result is even more 
surprising considering that the determination of Mx-my assumes 
some fixed parameters (e.g. the distance, the mass and the radius 
of the star, the abundance of the emitting plasma, and the temper- 
ature of the post-shock region) whose uncertainties are unknown 
and therefore not included in the errors on Mx-ray- Taking into 
account all these uncertainties, one would expect a scatter in the 
values of Mx-ny possibly larger than one order of magnitude. On 
the contrary, the evidence is that the observed scatter is within a 



factor » 2 and it can be explained only if the (unknown) un- 
certainties on the fixed parameters used in our calculations are 
smaller than expected. 

As for the almost constant value of Mx-my, we note that its 
estimate relies on the flux of the O vn resonance line and is based 
on the assumption of optically thin plasma. However it can be in- 
ferred that optical depth of the strongest emission lines produced 
by the plasma located in the post-shock region is non-negligible. 
In fact for a post shock region with dimension ~10 9 cm, filled 
with plasma at T=2 MK and n e =10 n cm 4 , the optical depth of 
the O vn resonance line is ~ 10. Supporting this scenario is some 
evidence of optical depth effects (i.e. the emitting plasma is not 
totally optically thin) which has been detected in the X-ray spec- 
trum of MP Mus (Argiroffi et al. 2009 ), implying that the O vn 
flux and the derived Mx-my could be underestimated. 
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The optical depth effects can be addressed by comparing 
mass accretion rates measured from the resonance line of the 
O vn triplet with the accretion rate measured using the flux of 
the intercombination plus the flux of the forbidden line of the 
same triplet, which are less affected by optical depth effects. We 
performed this comparison in the case of TW Hya and MP Mus. 
We found no significant differences in the former case, while the 
Mx-rm of MP Mus measured from the intercombination plus the 
forbidden line is twice as large as the Mx-my measured from the 
resonance line. The lack of high signal-to-noise spectra does not 
allow us to further investigate this issue using this approach for 
the other stars in our sample. However, this preliminary result 
suggests that the absorption from the stellar atmosphere is the 
main cause of the discrepancy between the optical and the X-ray 
mass accretion rates. 

Brickhouse et al. ( 120 1 Oi l analyzed X-ray spectroscopic data 
of TW Hya with high signal-to-noise ratio. This kind of analysis 
is limited to the case of TW Hya, which is the brightest CTTS 
in the X-ray and has been observed with Chandra for ~500 ks. 
These authors found that the flux and the densities measured 
from the Ne ix triplet agree with the accretion shock parameters 
derived from optical emission lines, while the density and flux 
derived from the O vn triplet disagree. They proposed therefore 
a model of "accretion-fed corona" in which the X-ray emission 
originates from three plasma components: a hot (T e « 10 MK) 
corona, a high density (« e « 6.0 x 10 12 cirT 3 , T e w 3.0 MK) 
post-shock region close to the shock front, and a cold less dense 
(n e ~ 2 x 10 11 cur 3 , r e w 2.0 MK) post-shock cooling region, 
with 300 times more volume and 30 times more mass than that 
of the post shock region itself. On the other hand, it is worth to 
note that, recently, Sacco et al. (1201 Oi l have shown that, if accre- 
tion streams are inhomogeneous and the chromospheric absorp- 
tion is efficient, their hydrodynamic model of accretion shock 
predicts that different He-like triplets measure different densi- 
ties of the X-ray emitting plasma. This is explained because the 
effects of absorption increases with wavelength and, as a conse- 
quence, the Neix (13.45 A) emission results to be less absorbed 
than O vn (21.60 A) emission. In particular, Sacco et al. (2010) 
showed that, in an inhomogeneous stream affected by chromo- 
spheric absorption, the average density measured with the Ne ix 
triplet is, in general, larger than that measured with O vn triplet 
although Ne ix lines form at temperatures higher than O vn lines 
(as noted by Brickhouse et al. I20101 >. Unfortunately, due to the 
lack of high signal-to-noise data covering the Ne ix triplet, we 
are not able to test this model for the other stars of our sam- 
ple. Further long exposure X-ray spectroscopic observations of 
the stars in our sample are required to test this model as well as 
optical depth effects on the estimate of M X -ray 

If the results presented in this paper are confirmed, we sug- 
gest here that the optical depth effects in the O vn resonance line 
may explain the almost constant Mx-my derived in our sample if 
the effect is larger in stars with larger optical derived mass accre- 
tion rates. Finally, it is worth to note that our analysis does not 
allow us to conclude that there is no relation between (Mopt) and 
Mx-my From our analysis we can just conclude that the relation, 
if exists, is masked by the (small) scatter in the values of Mx-ray 
To this respect, it is interesting to note that the lower right panel 
in Fig. |4] seems to suggest that Mx-my slightly increases for in- 
creasing (Mopt) if we consider that the value of Mx-my for Hen 
3-600 shows the largest error in X-rays. 



5. Conclusions 

We have carried out the first homogeneous comparison of optical 
and X-ray derived mass accretion rates for a sample of CTTSs. 
We compare the different optical tracers to each other, and to the 
X-ray derived accretion rates. We also analysed the variability of 
the CTTS TW Hya. Our findings lead to the following conclu- 
sions: 

1 . The mass accretion rates derived from the different optical 
tracers agree within the errors for each source, albeit with a 
large spread of typically « 1 order of magnitude (see Fig.[T|i. 

2. Ha 10% full width, whilst known to be a good indicator of 
accretion, is not good at measuring mass accretion rates. 

3. For the CTTS TW Hya (the only source within the sample 
for which a time analysis was possible), we find little vari- 
ation in the mass accretion rates (within each emission line 
tracer) over a period of 5 months (see Fig. |3j. The variabil- 
ity is much smaller than the range of accretion rates derived 
from different accretion tracers. 

4. The X-ray mass accretion rates are always smaller than the 
optically derived mass accretion rates for all sources within 
the sample (see Fig. |4). This can be explained if the ac- 
cretion streams are inhomogeneous and/or multiple streams 
with different densities are present simultaneously. In these 
cases, Sacco et al. (1201 Oi l have shown that the chromospheric 
absorption triggers a selection of the X-ray emitting shocks, 
absorbing preferentially the X-ray emission from high den- 
sity plasma components. As a result, only the low density 
plasma component of the post-shock plasma can be observed 
in the X-ray band, leading to a systematic underestimate of 
the mass accretion rate. 

5. We find that the X-ray derived mass accretion rate ranges 
within a factor of « 2 around 2 x 10~ 10 M yr~' (see Fig.|4]i, 
despite the fact that the range of optical mass accretion rates 
span almost 3 orders of magnitude. Some evidence of non- 
negligible optical depth of emission lines produced by post- 
shock accreting plasma (e.g. Argiroffi et al. 2009) may ex- 
plain the almost constant Mx-my if the effect is larger in stars 
with larger optical mass accretion rates. This issue deserves 
further investigation in future studies to assess the severity of 
optical depth effects on the estimate of mass accretion rates 
in the X-ray band. 
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